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Abstract 



A Darboux-transformed surface gravitational acceleration of the constant 
gravitational acceleration for a body endowed with an atmospheric layer is 
shown to turn the atmospheric free fall with quadratic resistance in the oppo- 
site motion, i.e., a free rising. Although the atmosphere of such a body may 
look completely normal, it is the time dependence of its gravitational field 
that produces this type of motion. This result is a consequence of general, 
one-parameter-dependent Darboux transformations in mathematical physics. 



I. INTRODUCTION 

The following time- dependent gravitational acceleration 



9{t;^) =9 



2— ln(/oi(t) + A) 



(1) 



where /oi(t) = Jq cosli^xdx and A is a parameter of the gravitational force, turns falling 
through a quadratic resistive - earth-like - atmosphere into just the opposite motion, i.e., a 
rising that for t — > cxd is of constant velocity. Eq. (1) is in fact the one-parameter Darboux- 
transformed acceleration of the constant acceleration g at the surface of a usual planet 
such as Earth. The result is a consequence of a mathematical scheme that has been called 
strictly isospectral Darboux technique (SIDT) and has been applied extensively by one of 
the authors to various fields of physics. Briefly, SIDT is a three step procedure that for 
the resistive atmospheric free fall means the following, 
(i) One starts with the equation 

representing the normal resistive falling motion through an atmosphere of constant friction 
coefficient per unit of mass e. 
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(ii) One shifts to a motion of the form 



'^'^ g2{t)-ev\ (3) 



dt 
where 

(72(t) = (?[-! + 2 tanh2(yi^t)] . (4) 

One interpretation of Eq. (3) is that the motion is upwards through a fluid medium 
producing a driving force proportional to the square of the velocity of the particle (one may 
call it antifriction) plus another force g2{t) of driving character for t < T = -^Arctanh-^ 
and turning friction-like afterwards. 

The first two steps are connected to each other through the particular solution of the 
initial normal motion as given by Eq. (7) below. 

(iii) The third step is the return back to an equation similar to the initial one, namely 

I = 9{t; A) - ev^ , (5) 

where the time-dependent forces g{t; X) are given by Eq. (1). The second and the third 
steps are also connected to each other through the general solution of the second problem 
(see Eq. (16) below). Thus, although the motion given by Eq. (5) looks normal, its solution 
v{t) corresponds to a rising on the vertical rather than a fall. Exactly as in the falling case, 
the rising has a short accelerating stage followed by a uniform rising with constant limiting 
velocity. Indeed, for sufficiently large times, the driving force is at its maximum value —eV"^, 
whereas g2{t oo) = +g is a constant friction force and one gets the value of V from the 
balance of the two forces. 

Here, we first review the usual resistive atmospheric falling according to a 40-year old 
discussion of Davis. Next, we briefiy describe the motion in a Darboux-transformed 
gravitational field. Some concluding remarks end up the work. 



II. DAVIS' RESISTIVE FREE FALL 



The free fall in quadratic resistive media has been studied by Davies in his book as a simple 
application of the Riccati equation for the falling velocity. Davis wrote the "readily found" 
solution of Eq. (2) in the following form 



vit) = v r'^''''y^' ] , (6) 

^ ' \1 + uotanh^/agt I ^ ^ 



where uq = vq/V and V = yg/e is the so-called limiting (terminal) velocity. Notice that 
Eq. (6) can also be written as 

v{t) = ytanh(ye^t + (3) , (7) 
where the arbitrary phase f3 is fixed through the initial condition, (3 = Arctanhwo. 
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Davis obtained the faUing distance sd(^) by integrating the velocity: 

soit) ^ [' v{t)dt^— In (cosh ^ + '^ sinh ^) . (8) 
Jo g \ V V V / 

The comparison with a set of 25 experimental data of parachute falls in the atmosphere of 
the earth was presented by Davis for the particular case uo = 0, i.e. Vp — Ftanh ^/egt and 

V"^ / at\ 

sd = — In cosh— ^Vt ln2 , (9) 

g \ VJ g 

where the last approximation is for large times. The motion sd is asymptotically a uniform 
falling with the limiting velocity V. 



III. ONE-PARAMETER DARBOUX TRANSFORMED FORCES AND 

RESULTING MOTION 



We now provide details of the Darboux constructions. We recall that nonrelativistic su- 
persymmetric quantum mechanics is a simple application of Darboux transformations, with 
a huge publication output during the last two decades. One usually starts with a Riccati 
equation, such as Eq. (2), that we write in dimensionless form 

n,: ^ + u'^l, (10) 
or 

having ui = tanh(r + 7) as solution. Next, it is easy to see that the right hand side 
corresponds up to a sign to a constant potential function for the Schrodinger equation at 
zero energy 

The particular solution wi = cosh(r + 7) is usually called a zero mode. It is connected to 
the Riccati solution through Ui = -^dwi/dr. Next, changing the sign of the first derivative 
in IZi one calculates the outcome gfij) using the same Riccati solution -Ui, i.e., — ^ + u\ = 
2tanh^(r + 7) — 1 := gfir). Thus, one considers a second Riccati equation 

112: — - = 2tanh^(T + 7) - 1 , (12) 
or 

of which we already know that it has ui as solution. The corresponding zero mode W2 — 
1/ cosh(r + 7) fulfills 

( d \ ( d \ d^W2 
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We are now interested in the general solution Ug2 of 1^2- To find it, one employs the 
Bernoulli ansatz Ug2{T) = «i(t) — j^, where /(r) is an unknown function. One obtains for 
the function /(r) the following Bernoulli equation 

^ + 2f{T)m{r) = l, (14) 

that has the solution 

/M^^. (15) 

where Xoi(r) = w\{y)dy, and we consider A as a positive integration constant that is 
employed free parameter of the force field. 

Thus, the general Riccati solution of TZ2 is a two-parameter function Ug2{T; 7, A) of the 
following form 



d 



d r f cosh(r + 7) 



dr 



Mr;/3;A)^- In [ ^ ^ J = :JZ i T-T:^urVZ^rVT7^VT^ ) ' (16) 



isinh2(r + 7) + 5(r + 7) + A^ 



where in the last A we absorbed the lower limit of the integral. The range of the A parameter 
is conditioned by Xo(r) + A 7^ in order to avoid singularities. According to SIDT one can 
use Ug2 to calculate a family of force functions gx := + u^2- Thus 

^A = l-2^1n(Joi(r) + A) . (17) 

This expression coincides with Eq. (1) when passed to dimensional quantities. Thus, there 
is a third Riccati equation, namely 

^ dU r, 

which is similar in form to TZ\ but possessing the solution Ug2. The corresponding linear 
equation is 

( d \ ( d \ d^wz _ ^ 

+ U3 = ^ - ^.^3 = , (19) 



where 



In the limit A — > 00, Eq. (18) goes into Eq. (10) because Ug2 — > ui and g\ — > 1. Thus, one 
can say that Eq. (18) is the one-parameter SIDT generahzation of the Riccati equation (10). 

The problem of imposing initial conditions can be solved similarly to Davis. Since the 
solution Ug2 has two parameters, 7 and A, we shall fix the phase parameter through the 
initial condition Mg2(0) = uq and the other parameter. This leads to the following cubic 
algebraic equation 

tanh^ 7 — Mo tanh^ 7 — tanh7 = —{uq + —) . (21) 

A 
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The three solutions of this equation can be obtained by Cardano's formulas and any of them 
can be used to fix 7. The real one is the most simple 



1 / EV3_^-i/3\ 
71 = Arctanh- Iuq ^ 1 , (22) 

where E = ( ^+^^f +^' ) A = -3X^ - X^uq, and B = 27X^ + 27X^uo - 9X^uo - 2X^ul The 
other two solutions are complex and being more complicated will not be reproduced here. 
Let us calculate the distance in the motion 7^3. It is given by 

S3(r) = r u,2 = In f ^7l'"l^ l ^ 1^ Ue^e-'A ^-Vt+^ln2 (23) 
Jo \Joi{t) + XJ ^ ^ g 

that shows that in the large-time asymptotics the motion is a vertical rising of constant 
velocity V . 



IV. CONCLUSION 



One may say that if there exist planetary-like astrophysical bodies with gravitational fields 
of the type given by Eq. (1), then an anti-parachute effect might occur in their atmospheric 
layer. It is hard to believe that nature can produce such fields. For example, in the case of 
the common SIDT application to supersymmetric quantum mechanics, the scheme is used 
merely to extend the class of exactly solvable potentials and in general is considered to be a 
curiosity. However, we point out the debate on the modifications of the gravitational fields 
in superconductors as a consequence of macroscopic quantum effects with applications to 
neutron stars. |0] 

Finally, we would like to make the following comment. SIDT introduces singularities in all 
the transformed quantities. The parameter A gives the measure of the time interval where 
the singularity is located in the past (t < 0). A small A means a singularity close to the 
moment of parachute jumps and the time dependence is significant. On the other hand, a 
large A leads to very small variation in time of the planetary gravitational field that can be 
thought of as a cosmological evolution implying a singularity very far in the past. 



The authors are much indebted to Dr. M. Nowakowski for discussions and active criticism. 
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Appendix 

In general, for a Riccati equation of the type 

dy 

— = ay^ + , a ,(3 = const. (Al) 
dx 

the solutions can be written as follows 
(i) for I > 

y = ]j~^tan (J^px + (f)i\ or y = cot (-J^fSx + cj)'^ {A2) 



[ii) for I < 



y = ]J -^tanh (^^j^px + ipij , or y = ]j -(^coth (^^j^px + ip2j . {A3) 

Comparing the Newtonian falling case with Eq. (Al) one gets 

a — e 

a = -e, p = g =^ =_<0, (AA) 

^ 9 

and therefore the solution is given by 

v(t) = J- tanh {Jtgt + iff) , (A5) 
V e 

which is Davis' result. The falling velocity v{t) presents a substantial time dependence - 
acceleration stage - only for a short period and rapidly turns into the limiting constant 
velocity V — '- 
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3: 



The rising distance S3 for the same 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
tau 



Fig. 5 The zero modes wi, W2, and ws{X = 




Fig. 6 The normalized force g\ for 71 = 1 
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